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A high-precision numerical calculation is reported for the self-energy correction to the hyperfine 
splitting and to the bound-electron g factor in hydrogenlike ions with low nuclear charge numbers. 
The binding nuclear Coulomb field is treated to all orders, and the nonperturbative remainder 
beyond the known Za-expansion coefficients is determined. For the 3 He + ion, the nonperturba- 
tive remainder yields a contribution of —450 Hz to the normalized difference of the IS and 2S 
hyperfine-structure intervals, to be compared with the experimental uncertainty of 7f Hz and with 
the theoretical error of 50 Hz due to other contributions. In the case of the g factor, the calculation 
provides the most stringent test of equivalence of the perturbative and nonperturbative approaches 
reported so far in the bound-state QED calculations. 

PACS numbers: 31.30.jf, 31.30.js, 32.10.Fn 
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The hyperfine structure (hfs) of the ground state of 
hydrogen is experimentally known with a relative accu- 
racy oflxl0~ 12 [lj, this measurement having for a long 
time been among the most precise ones in physics. One 
of the remarkable features of the hfs is an important role 
of the binding effects in its theoretical description. For 
the self-energy (SE) correction to the hfs, the binding 
effects change the sign of the correction already for a 
nuclear charge number Z = 8 and make the expansion 
in the binding-strength parameter Za completely inade- 
quate for high values of Z (here, a is the fine-structure 
constant). Large coefficients of the Za expansion and the 
high accuracy of experimental results call for an all-order 
(in Za) approach in the theoretical description of the hfs 
even for systems as light as hydrogen. 

The high-precision all-order calculation of radiative 
corrections for hydrogen is a notoriously difficult prob- 
lem. This point can be illustrated by considering the SE 
correction to the Lamb shift. Its accurate evaluation to 
all orders in Za was first accomplished by P. J. Mohr in 
1974 0] for Z > 10, while an analogous calculation for 
Z = 1 was not realized until two decades later Q. 

All-order calculations of the SE correction to the hfs 
started in late 1990s [J, [a @- The first attempt at a 
numerical evaluation for hydrogen was made at the same 
time in Ref. [7|. Due to insufficient numerical accuracy 
at Z = 1 in that work, the goal was reached in an in- 
direct way: the known terms of the Za expansion were 
subtracted from the all-order numerical results for Z > 5, 
and the higher-order remainder thus inferred was extrap- 
olated towards Z = 1. The result obtained was used as 
an important theoretical input for the determination of 
the muon mass from the muonium hfs measurements |8| . 

The accuracy of the numerical evaluation of the SE 
correction to the hfs was improved by several orders of 



magnitude during the past years [9|, LLIJ- However, the 
precision obtained was still insufficient for a direct deter- 
mination of the higher-order SE remainder at Z = 1, 
and an extrapolation procedure was again employed. 
The studies [9j, [l0( provided a remainder value for the 
normalized difference of the IS and 25* hfs intervals, 
A 2 = 8AE 2 s — AEis [IH, in 3 He + and demonstrated a 
2<7 deviation of the theoretical prediction from the exper- 
imental result [l^LUI- The accuracy of the extrapolation 
procedure of Refs. [3, [l0| has recently become a sub ject 
of some concern. In particular, there is an opinion [14j | 
that the error of the extrapolation is four times larger 
than given in Refs. [9j, llOj , which would bring theory and 
experiment back into agreement. 

The main goal of the present investigation is to perform 
the first direct, high-precision theoretical determination 
of the higher-order remainder of the SE correction to the 
hfs of the 15 and 2S states of light hydrogenlike ions. In 
addition, we carry out a related study of the SE correc- 
tion in the presence of an external homogeneous magnetic 
field, i.e., the SE correction to the g factor of an electron 
bound by a spinless nucleus. 

High-precision experimental investigations for the 
bound-electron g factor have a shorter history than those 
for the hfs but are not less important. A relative accuracy 
of 5 x 10~ 10 was reached in recent microwave measure- 
ments in hydrogenlike carbon and oxygen [15|, [l6| , thus 
provi ding a new tool for the determination of the electron 
mass [ItJ ■ A recent proposal [18[ to employ laser spectro- 
scopic techniques in these measurements opens perspec- 
tives for improving the experimental accuracy (particu- 
larly, for the helium ion) up to the level of 10~ 12 . 

Already at the present level of experimental accuracy, 
the theoretical description of the bound-electron g factor 
has to be performed to all orders in Za. The numerical 
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FIG. 1: The electron self-energy in the presence of a magnetic 
field. The double line indicates the bound electron and the 
wave line with a cross is the magnetic field. 



precision often becomes a matter of crucial importance 
in such calculations. So, an increase of the numerical ac- 
curacy for the SE correction to the g factor by an order 
of magnitude achieved in Ref. [lj| as compared to the 
previous evaluations [y, |2(J, UM resulted in an improve- 
ment in the electron mass value. In order to match the 
10~ 12 level of accuracy anticipated in future experiments 
on the helium ion, the precision of numerical calculations 
of the SE correction should be enhanced by several or- 
ders of magnitude. This task will be accomplished in the 
present work. 

The SE correction in the presence of an external (mag- 
netic) potential V m:ign is graphically represented by two 
topologically nonequivalent diagrams in Fig. [TJ Formal 
expressions for them can be obtained by considering a 
first-order perturbation of the SE correction to the Lamb 
shift by Knagn- Perturbations of the reference-state wave 
function, the binding energy, and the electron propaga- 
tor give rise to the irreducible, the reducible, and the 
vertex contributions, respectively. General formulas for 
these contributions are known and can be found in our 
previous study 1 10] ; for a detailed analysis we direct the 
reader to Ref. [22|. The irreducible part reads 



AE ir = (a\Z(e a )\6a) + (8a\Z(e a )\a) , 



(1) 



where S(e a ) is the SE operator defined so that its diag- 
onal matrix element (a | £(£„)! a) yields the one-loop SE 
correction to the Lamb shift [2J, and \Sa) is the first-order 
perturbation of the reference-state wave function \a) by 
Knagn- The reducible part is given by 

d 

AE Ted = {a\V magn \a) (o| — E(e)|o> , (2) 

and the vertex part is 

■ 2 poo 

AE veI = — / dco Yl 

^ j-oo „^r 2 

[e a —v -ue ni )(e a - uj - ue n2 ) 

where a M = (l,ot) are the Dirac matrices, u = 1 — iO, 
and D^ v is the photon propagator. 

The calculation of the irreducible part is similar to the 
evaluation of the diagonal matrix element of the SE op- 
erator. It is performed here by a generalization of the 



approach of Ref. 3] , with the use of the closed- form an- 
alytical representation for the perturbed wave function 
| Sa) [23|. The evaluation of the reducible and vertex 
parts is more difficult. It is carried out after splitting 
them into several parts, 



AE ied = AE, 
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(4) 
AEiH + AEiH + A2S&+> , (5) 



where the upper index (a) labels the contributions in- 
duced by the reference-state part of the electron propa- 
gators and the other indices specify the total number of 
interactions with the binding field in the electron prop- 
agators [the index (i+) labels the terms generated by 
>i such interactions]. The reference-state contributions 

A£^ed an d AEvtr are separately infrared divergent. The 
divergences disappear when the contributions are regu- 
larized in the same way and evaluated together. The 
zero-potential parts AE^l and AE^d are separately ul- 
traviolet divergent. They are covariantly regularized by 
working in an extended number of dimensions and cal- 
culated in momentum space. The remainder of the re- 
ducible part A£ rcd contains at least one interaction 
with the binding field in the electron propagators and is 
finite. In its evaluation, advantage was taken of a general- 
ization of the numerical procedures originally developed 
in Ref. Q. 

The remaining vertex contributions AEvd and AEver 
contain three electron propagators and represent the 
most difficult part of the calculation. The key to the 
success was to isolate the one-potential vertex contribu- 
tion AE'vor an d to calculate it without any partial-wave 
expansion in momentum space. For the SE correction to 
the g factor, such a calculation has been carried out in 
Ref. (19j |. employing the fortunate fact that in momen- 
tum space, the interaction with the homogeneous mag- 
netic field is expressed in terms of the (gradient of the) 
(5-function. This is not the case for the hfs, and the cal- 
culation of this contribution is much more difficult. The 
general expression for A_Evor reads 
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dp dp' dp" 



d A k V c {p") 



(2tt) 9 J (2tt) 4 k 2 

x V>«(P) l<yS{P ~ k ) 7o5(p - p" - k) 

x loV magn (p'") S(p' - k) 7 - MP') , (6) 

where p"' = p - p' - p", p a = p' = e a , pi = p{," = 0, 
S(p) = l/(j/) — m) is the free electron propagator, and Vc 

is the Coulomb potential. Effectively, AEvd is a two-loop 
contribution because two momentum integrations (over 
d 4 k and dp") need to be performed analytically. They 
are carried out after joining denominators by introduc- 
ing 4 auxiliary Feynman parameters. Next, we integrate 
over all angular variables except p p' , which leaves 7 in- 
tegrations (3 over the kinematic variables and 4 over the 
auxiliary parameters) to be carried out numerically. The 
numerical evaluation is rather time-consuming (about a 



TABLE I: SE correction to the hfs of nS states of hydrogenlike 
ions. SE n = AE n /[(a/-K) E F (n)] and F n is defined in Eq. Q. 



8E n 



F„ 



Ref. 



0.438 fOl 842 (2) 



0.373 467 603(3) 
0.307 583 837(4) 
0.241 005 729 (6) 
0.174 026 210(7) 



-13.8308(42) 

-13.8(3) 

-15.9(1.6) 

-12 (2) 

-14.1159(10) 

-14.4120(4) 

-14.6962(2) 

-14.9673(2) 



[10] 

[25] 



0.438 692 275 (3) 



0.375 352 040(4) 



0.311203192(5) 
0.246 665 422(7) 
0.181938 683(10) 



-6.1205(84) 
-6.2 (9) 
-7.8(1.4) 
-6.9129(11) 
-6.9 (4) 
-8.2 (9) 
-7.5833 (5) 
-8.1698 (3) 
-8.7069 (2) 



[10] 
[9] 

[10] 
[9] 



month of processor time for each value of Z) but the 
crucial point is that it does not involve any partial wave 
summations. 

(24-) 

The remaining vertex contribution A_Ever contains 
bound electron propagators, and so the partial-wave ex- 
pansion in its evaluation is unavoidable. However, the 
convergence of this expansion turns out to be very good, 
provided that the integrations over all radial variables 
are first carried out. For instance, at Z — 1, the sum 
of only the first two partial waves for the hfs yields a 
result with a relative accuracy of 10~ 5 . The good con- 
vergence is due to the separation of the one-potential 
contribution AEiJ introduced in this work. About 120 
partial waves included in the actual calculation and the 
extended-precision arithmetics employed allowed us to 
control the calculation to a level of 10 -9 . 

The results for the SE correction to the hfs can be 
conveniently represented as 



AE„ 



E F (n) — a 00 + (Za) a w + (Za) 2 (L 2 a 22 + La 2 



+ a 20 ) + (Za) 3 La 31 + (Za) 3 F n (Za) 



(7) 



where Ep(n) is the non-relativistic hfs value, L = 
ln[(Za) ], and o^ are coefficients of the Zcv-expansion 
known today: a 00 (nS) = 1/2, a w (nS) = -8.03259003, 
a 22 (nS) = -2/3, a 2 i(lS) = -1.334504, a 21 (2S) = 
0.317104, 020(15*) = 17.122339, a 20 (2S)_ = 11.901105, 
a 31 (nS) = -13.307416, see recent works [Hill] and ref- 
erences therein for earlier studies. F n is the higher-order 
remainder, which should addressed in a numerical all- 
order approach. 

The results of our numerical calculation of the SE cor- 
rection to the hfs of the IS and 25* states of light hy- 
drogenlike ions with Z < 5 are presented in Table [U 



TABLE II: SE correction to the IS bound-electron g factor, in 
units of 10 -6 (ppm). Hi is the higher-order remainder defined 
by Eq. ((SJ and obtained by a direct evaluation; Hijextr.) 
denotes the extrapolated results. The results of Refs. [191. [211] 
are scaled for the present value of a [g]. 



Affi 



Hi 



Hi (extr. 



Ref. 



2 322.840 245(1) 
2 322.840 3(1) 
2 322.840 2(9) 
2 322.904052(4) 
2 323.014310(8) 
2 323.175 54(2) 
2 323.392 99(2) 



12(31) 



23.1(2.8) 
22.88(70) 
22.57(30) 
22.35(16) 



23.39 (80) 



23.03 (44) 



[19] 

[21] 



The fine-structure constant of a" 1 = 137.03599911 [8j 
was employed in the calculation. Since the current un- 
certainty of a (3 ppb) does not influence the numerical 
accuracy of the higher-order remainders, a is assumed 
to attain exactly the value indicated. Good agreement 
is observed with the extrapolated values of the higher- 
order remainder obtained previously [7|, [9j, [l0| and with 
the Za-expansion result of Ref. [25j , but their accuracy 
is increased by several orders of magnitude. 

Our calculation removes a significant source of uncer- 
tainty in the theoretical predictions for the normalized 
difference of the IS* and 2S hfs intervals in hydrogen and 
helium-3 ion, A 2 = 8AE 2 - AE 1 . For 3 He+, the SE re- 
mainder determined in this work amounts to —0.450 kHz. 
Combining this result with other theoretical contribu- 
tions to A 2 described in detail in Refs. [Tlj, [2|| |27I ]. 
we obtain the total theoretical value A 2 ( 3 He + ) t hco = 
— 1190.135 (50) kHz, to be compared with the experimen- 
tal result A 2 ( 3 He+) cxp = -1189.979(71) kHz [H, [l|. 
Our calculation of the SE remainder improves the accu- 
racy of the theoretical prediction by a factor of three, 
as compared with Ref. [2g]. For hydrogen, the the- 
oretical and experimental results read A 2 (H) t heo = 
48.9541(23) kHz and A 2 (H) cxp = 49.13 (13) kHz 0,H|, 
correspondingly. 

For the g factor, the results of our numerical evaluation 
can be parameterized as 



Ag n = - 1 + (Za) 2 b 20 + (Za) 4 (L 641 + b. 



40) 



+ (Za) 5 H n (Za) , (8) 



where bij are known coefficients of the Za expan- 



\n- 2 
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b\\(nS) 



.VI 



n~ s ) _bA {lS) 



sion: b 2 o(nS) 

-10.236524, 640(25) = -1.338464, see Ref. [29( and ref- 
erences therein. H n is the remainder incorporating all 
higher-order contributions. It is remarkable that for the 
g factor, the higher-order remainder enters in the relative 
order (Za) 5 rather than in the relative order (Za) 3 , as in 
the case of the hfs. This means that cancellations in ex- 
tracting the remainder from numerical results for Z = 1 
are by four orders of magnitude larger for the g factor 
than for the hfs. 



The results of our numerical calculation of the SE cor- 
rection for the 15 bound-electron g factor are presented 
in Table [ill For hydrogen, they are consistent with val- 
ues reported previously pjj, [21j but are by two orders 
of magnitude more accurate. At the same time, all ten 
digits of our numerical all-order result for Z — 1 coincide 
with the value obtained within the Za-expansion. The 
fact of this coincidence can be considered as one of the 
most stringent tests of consistency of the two main the- 
oretical approaches presently developed in bound-state 
QED. Similar agreement between the Za-expansion and 
the all-order approach was observed for the 25 state; the 
corresponding results will be presented elsewhere. 

The accuracy of the direct numerical determination 
of the g-factor remainder Hi for Z = 1 and 2 can be 
increased by extrapolating our results obtained for the 
higher-Z region. We employ the extrapolation proce- 
dure described in Refs. [iJLLll an d the numerical data for 
the remainder H\ for Z as high as 20 in order to obtain 



the improved results listed in Table [TT| under the label 
ff^extr.). 

To conclude, we have performed high-precision all- 
order calculations for the SE correction to the hfs and 
to the bound-electron g factor in light hydrogenlike sys- 
tems, improving the numerical accuracy by several orders 
of magnitude as compared to the previous evaluations. 
Accurate nonperturbative results have been obtained for 
the higher-order SE remainder for the hfs. We remove 
an important source of uncertainty in theoretical predic- 
tions for the normalized difference of the 15 and 25 hfs 
intervals in hydrogen and the hclium-3 ion and increase 
the theoretical accuracy by a factor of three. 
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